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Abstract 

We present and analyze a discontinuous variant of the Zip- version of the boundary clement 
Galerkin method with quasi-uniform meshes. The model problem is that of the hypersingular 
integral operator on an (open or closed) polyhedral surface. We prove a quasi-optimal error 
estimate and conclude convergence orders which are quasi-optimal for the /i-version with 
arbitrary degree and almost quasi-optimal for the p- version. Numerical results underline the 
theory. 

Key words: Zip-version with quasi-uniform meshes, boundary element method, discontinuous 
Galerkin method, hypersingular operators 
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1 Introduction 

Discontinuous approximations of solutions to boundary value problems of second order are well 
established and have their advantages in comparison to conforming methods, e.g., by providing 
more flexibility in the design of discrete spaces. This applies not only to finite elements but 
also to the boundary element Galerkin method (BEM), that is, finite element approximations 
of solutions to boundary integral operators. Whereas there is a huge amount of literature about 
discontinuous finite elements, relatively little is known on non-conforming approximations of 
boundary integral operators. In this paper we analyze a discontinuous Zip-approximation of 
the hypersingular operator governing the Laplacian. In principle, our analysis applies to other 
hypersingular operators as well, e.g. from linear elasticity and acoustics (Helmholtz equation), 
in the sense that the main tools from fractional order Sobolev spaces remain identical. But some 
specific estimates are indeed non-trivial for these applications. Note also that boundary element 
approximations of solutions to problems governed by weakly singular operators or by equations 
of the second kind can use discontinuous approximations, they are conforming. Hypersingular 
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operators are the only boundary integral operators from second order elliptic problems that are 
challenging when one tries to relax the continuity of basis functions. 

We started our study of non-conforming approximations of hypersingular operators with the 
analysis of Lagrangian multipliers |9j (for the implementation of essential boundary conditions) 
and Crouzeix-Raviart elements [13] . In [10] we extended the Lagrangian multiplier technique to 
a domain decomposition method, and in [5] we presented a domain decomposition variant based 
on the Nitsche coupling that avoids the additional space needed for a Lagrangian multiplier. All 
these results are only about lowest order approximations. Indeed, their analysis makes heavy 
use of finite-dimension arguments (like norm equivalences and inverse properties) and Sobolev 
norm estimates on fixed domains (sub-surfaces), and do not extend to high order methods. 

In this paper we study the Zip- version with quasi- uniform meshes. Key point of our analysis 
is to avoid finite-dimension techniques from previous papers and to work with general Sobolev 
norm estimates (not restricted to finite-dimensional spaces) so that on elements there is no 
restriction for polynomial degrees. In this way, we are able to prove quasi-optimal convergence 
of a discontinuous boundary element Galerkin method (Theorem 3.1) that holds on conforming 
quasi- uniform meshes for general non- uniform polynomial degree distributions. The small price 
to pay is that this Cea estimate involves three different Sobolev norms (L 2 , H 1 / 2 and H s 
with s > 1/2). Nevertheless, it leads to a quasi-optimal estimate in the case of the h- version 



(Corollary |3.2[ ) and an almost quasi-optimal estimate for the p- version (Corollary 3.3). In the 
latter case, almost quasi-optimal means that there is a log 3//2 (p)-perturbation of the quasi- 
optimal estimate for the conforming method. Though, such logarithmic perturbations are known 
from the previous results (except for Crouzeix-Raviart elements) and it is remarkable that our 
/i-estimate is free from this. 

Principal techniques used in the analysis are of the Strang type (discrete ellipticity and 
consistency of the discrete formulation). In this sense, the approach follows the same lines as 
our analysis of a Nitsche domain decomposition |5j, nevertheless relying on local instead of global 
estimates. For ease of presentation, we assume that meshes are conforming. However, techniques 
from domain decomposition in [5] can be used without difficulty to extend our analysis to meshes 
which are conforming only on sub-domains (sub-surfaces). In this case, an additional logarithmic 
perturbation appears in the error estimate. 

The remainder of this paper is structured as follows. In the next section we define some 
Sobolev norms, present our model problem (on a closed polyhedral surface), specify the standard 
boundary element (Galerkin) method, and recall an integration-by-parts formula for the hyper- 
singular operator. In Section [3] we present our discontinuous boundary element method, state 
the main results (quasi-optimal convergence, Theorem |3.1|) and conclude convergence orders for 



the lowest-order h- version (Corollary 3.2), the p- version (Corollary 3.3), and the h- version with 
arbitrary polynomial degree (Corollary 3.5). For presentation of the last two corollaries we need 
to recall results on the regularity of the solution to our model problem (in terms of appearing 
singularities) which is done in the same section. Technical details and the proof of the main 
theorem are given in Section [4j In Subsection 4.2 we discuss the changes which are necessary to 
analyze the model problem on an open polyhedral surface. In fact, there is some improvement 
for a single smooth surface piece. In Section [5] we present several numerical results that confirm 
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our error estimates. 

Throughout the paper, a < b means that a < cb with a generic constant c > that is 
independent of involved parameters like h or p. Similarly, the notation a > b and a ~ b is used. 



2 Sobolev spaces 

We consider standard Sobolev spaces where the following norms are used: For Q C M n and 
< s < 1 we define 

II ||2 _ || ||2 , | |2 

\\ u \\h 3 (q.) ■— H m Hl 2 (q) l u lfl" s (n) 

with semi-norm 

/ /" /" -u(y)| 2 , ,\V2 

For a Lipschitz domain f2 and < s < 1, the space H s (0,) is defined as the completion of Cq°(£1) 
under the norm 

/ , .2 f u (%) 2 , \ 1/2 

\U s„ m := \\U L sm + / tt: ; „_ dx 



For s E (0,1/2), || • [|jj s /Q\ and || • ||fl-s(n) are equivalent norms whereas for s G (1/2,1) there 

holds H S (Q) = Hq(Q), the latter space being the completion of C^°(0) with norm in H s (£l). 
For s > the spaces H~ S (Q) and H~ S (Q) are the dual spaces (with L 2 (f2) as pivot space) of 
H s (£l) and H s (£l), respectively. For norms of vector valued functions we use the same notation 
as for scalar functions. 

Let r be a piecewise plane, open or closed, polyhedral surface with faces Tj, j = 1,..,,L. 
Throughout the paper we will identify faces with polygonal subsets of IR 2 . First we will analyze 
the case of a closed surface, and in Section 4.2 we will mention particular changes which are 
necessary to analyze the case of an open surface. 

Our model problem is: For a given function f G L 2 (T) find u G H l / 2 (T) such that (u, l)r = 
and 

Wu{x) := -^j^t-t I u(y) „ 9 1 dS(y) = f(x), xeT. (2.1) 
4tt on(x) Jy on(y) \x — y\ 

Here, n is the exterior normal unit vector on T, and (•, -)r denotes the L 2 (r)-inner product and its 
extension by duality. Throughout, this generic notation will be used for other domains/surfaces 
as well, indicated by the index. 

A variational formulation of (2.1) is: Find u G H l / 2 (T) such that (u, l)r = and 

{Wu,v) T = (f,v) T W G H l l 2 (T). (2.2) 



A standard (i.e. conforming) boundary element method for the approximate solution of (2.2) 
is to select a piecewise polynomial subspace Hh p ,coni C H l l 2 (T) and 

H h P ,coni '■= i v e H hpiCon{ ; (v, l) r = 0}, 
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and to define an approximant u^cont £ H® conf by 

(Wu hpyCon{ ,v) r = (f,v)r Vv € H% pcon{ . (2.3) 

The bilinear form with hypersingular operator W is usually calculated by making use of its 
relation to a bilinear form with weakly singular operator V defined by 

V<p{x) :=]- [ dS(y), G (^-^(r)) 3 , x G T. 

4vr y r |x - y\ 

There holds the relation between the operators W and V: 

W = curlr (V curl r ) (2.4) 

as a linear continuous mapping from H l l 2 (T) to H^ 1 / 2 ^) (see [E]), so that 

{Wu,v)r = (curlr «, V curl r v)r Vu,i> G F 1/2 (r), 

cf. also j9j Lemma 2.3]. Here, curlr is the surface curl operator and curlr its adjoint operator. 
In the following, curlQ will denote the restriction of curlr onto Q C T and we will use curl^ 

as the surface curl (rotated gradient) on a reference element Q in local coordinates. 

Let us recall from [9j that integration by parts for curlr on any face Tj gives rise to a linear 
bounded operator 



(2.5) 



{v G H 1 ' 2 ^)- Wv G H-^iTj)} -> H-V 2 {8£j) 

v i-)- tj ■ (V curl r v)\ dTj 

with 

(Wv,w)r j = (curlrj w,V curlr u)r,- + (tj • V curlr v,w)gr j Vw G // 1 (r :) '). 

Here, tj is the unit tangential vector along the boundary of Tj, with mathematically positive 
orientation which is compatible with the direction of the normal vector n. In fact, it is easy 
to see that the operator (2.5) maps to H~ e (dTj) for any e > so that the test functions uu in 
(2.7) can be less regular (we don't make use of this here). Below, we need an integration-by- 
parts formula on sub-domains (sub-surfaces) Q C T. In order to have a well-defined bilinear 
form (Wv,w)q for any sufficiently regular subset Q of T we assume that Wv G L 2 (T). Then 
analogous arguments prove that for any Lipschitz surface Q CT with unit tangential vector tQ 
(with mathematically positive orientation) there holds 

{v G H l / 2 (T)\ Wv G L 2 (r)} -> H-V 2 (dQ) (2 g) 

v H> t Q ■ (V curlr v) \ dQ ! "'" 

with 

(Wv, w)q = (curlQ w, V curlr v)q + (tQ ■ V curlr v, w)qq ^w^H 1 (Q). (2.7) 
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3 Discontinuous Galerkin method 

Let Qh be a conforming quasi-uniform mesh which is compatible with the faces of T and whose 
(closed) elements are shape- regular triangles or quadrilaterals Q. The maximum diameter of the 
elements is denoted by h. The set of edges of Qh is denoted by £h- We also need the skeleton 
of the mesh 

For sufficiently smooth functions v on V we define jumps as follows. For e G £ h with e = QiHQ2, 

Qi,Q2 e Q h , 

[v]\e:=(v\ Ql -v\ Q2 )\ e . (3.1) 

Here, the selection of the numbering Qi, Q2 is arbitrary but fixed. 

We also assign unit tangential vectors, on the boundary of Q G Qh- On edges, unit 
tangential vectors are inherited (using the association e^Qi previously mentioned): 

t e :=t Ql \ e e G Eh- (3.2) 

We define the broken surface curl operator for sufficiently smooth functions by 

{cwv\ h v)\ Q := curl Q (v| Q ) VQ G Q h . 

We use discontinuous hp-spaces with variable degree for discretization. For a triangle Q G Qh 
let Pp(Q) denote the space of polynomials of total degree p, and for quadrilaterals let P P {Q) be 
the space of (bi) linearly transformed polynomials of degree p in both coordinates on a reference 
square. Throughout, we will use the generic notation Q for a reference element (the unit square 
for quadrilaterals and a fixed triangle for triangular elements) . 

For given polynomial degrees p = p(Q) (independent for different elements Q) we introduce 
discrete spaces 

X hp := {v G L 2 (r); v\ Q G P p(Q) (Q) VQ G Q h ) 

and the subspaces 

X° hp := {v G X hp - (v, l)r = 0}. 

The discontinuous Galerkin boundary element method then reads: 
Find Uh p G X® p such that 

a h (u hp ,v) = (f,v) r G X® ip . (3.3) 

Here, 

a h (v,w) := (y curb, curler + (Tv, [w]) lh - ([v],Tw) lh + u([v], [w]) lh (3.4) 
where v > is a given parameter. Furthermore, the operator T is defined as 

(Tt;)| e :=(Fcurl h t;)|e-te (e G £h)- (3.5) 
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It goes without further discussion that the integral mean zero condition in X® p can be imple- 
mented by a rank one Lagrangian multiplier. 

Note that curlew G L 2 (T) for v G X^ ip so that the tangential trace of G H 1 ^) is 



well defined without employing formula (2.7). 



The analysis will be based on "broken" Sobolev (semi-) norms 

'^(s,) : = E M^(q), h\\ 2 H t {Q h ) ■= \ v \h( Qh) + ^WMWh^) + 1 / v \ 2 (1/2 <«<!)■ 



Theorem 3.1. Let u G H r (T) with r G (1/2,3/2) 6e i/ie solution of (2.2), and Zef > 0, 



(5 G (0,1/2). Then, the discrete problem (3.3) is uniquely solvable and there exists a constant 
C > 0, depending on r and 5, but not on u, u, the actual mesh and polynomial degrees, such 
that for any s G (l/2,min{r, 1 — 5}] there holds the quasi-optimal error estimate 

h-u hp \\ H i /2{Qh) <c 

inf I \\u — v\\ T7i/2ir\ H — T77T, , ;„ I h~ l l 2 \\u — v\\ + /i s_1//2 ||u — ull f/sfn ) >. 

u6 xo p nCO(r)l" / (r) vV 2 {a- 1/2) 3 / 2 V 11 " L (r) 11 MjH (i V J 

A proof of this result will be given at the end of Section [4j The parameter 5 above is needed 
to maintain a fixed positive distance of s to 1 so that norms involved in the proof are uniformly 
equivalent and the constant C does not depend on s. This will be needed for the estimate of 



the p- version when s — > 1/2, see Corollary 3.3 below. First let us consider the h- version with 



quasi-uniform meshes and piecewise (bi)linear functions. Standard approximation results prove 
that in this case the discontinuous Galerkin boundary element method performs as well as the 
standard (conforming) BEM. 

Corollary 3.2 (/i-version lowest degree). Let p = 1. For u G H r (T) with r < 3/2 there holds 

h ~ u hp\\ H y 2 (Q h ) ~ (1 + u ~ 1/2 ) hr ~ 1/2 \\ u \\Hr(r)- 

Fixing a mesh and improving the approximation by increasing polynomial degrees leads to 
the p-version. Convergence analysis for conforming methods and based on standard Sobolev 
regularity proves convergence orders which are analogous to those of the /i-version (regularity 
order minus 1/2), but without upper bound. In the case of polyhedral surfaces, regularity is 
limited and depends on the angles at edges and corners where singularities appear. Specific 
approximation analysis for these singularities show that convergence orders are twice that of the 
h- version when element boundaries are aligned with edges, as is the case in this paper. In order 
to claim and prove a precise error estimate we would need to recall the corresponding regularity 
results by Dauge, von Petersdorff and Stephan [3 [T71 [18] . In order to improve readability, and 
following the results from these references, we reduce this presentation to the assumption that 



the solution u of (2.2) can be written like 



u = u 



rcg 



7£{edges} ue{vertices} 7£{edges}, ue{vertices}, vey 
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where u 1 , u v and u yv are the edge-, vertex-, and edge- vertex singularities of u, and u reg is a 
remainder of higher regularity. The sets {edges} and {vertices} denote the edges and vertices 
of the polyhedral surface V. The singularities are of the types (ignoring cut-off functions and 
smoother parts) 

v? = |logdist(-, 7 )|'dist(-, 7 f (I G {0, . . . , s}, s G IN U {0}, fi > 1/2) 

u v = | log dist(-,v)|* dist(-,t;) A (t G {0, . . . , q}, q G IN U {0}, A > 0) 

u? v = | logdist(-, 7l )|' I logdist(-, 72 )|* dist(-, 7l ) A ^ dist(-, l2 y (I G {0, . . . , s}, t e {0, . . . , q}, 

H > 1/2, A > 0) 

where 7 i and 7 2 denote the edges that meet at the vertex v. The smallest parameters /i, A, and 
largest number s, q limit the convergence order of the p-version. 

The following result shows that the p- version of the discontinuous Galerkin BEM converges 
almost as fast as the conforming version: the error estimate behaves the same asymptotically 
except for an additional log 3 / 2 p factor which, however, can be compensated by choosing v ~ 
(logp) 3 . 

Corollary 3.3 (p- version) . Let Xh p be defined by a fixed mesh that consists of triangles and/or 
parallelograms and let p denote the minimum polynomial degree. We assume that the solution 
u of (2.2) is of the type (3.6) with minimum parameters [i and A and maximum parameters q, 
s and that the regular part u reg is as smooth as needed, which can be achieved by considering 
sufficiently many singularity terms. Then there holds 

\\U - U hp \\ Hl „ [Qh) < (1 + V- 1 ' 2 log 3 / 2 (l + p))p -2min { A + l/2, M } + p)> 

where 

s + q + 1/2 if X = fi- 1/2, 



P 



s + q otherwise. 



Proof. By Theorem 3.1 there holds for a fixed mesh 

\\ u - u h P \\ H i/2 (0h) < inf (\\u-v\\ m/2(r) +u~ 1/2 {s- l/2)" 1 ||n-i;||^ (r) 

and by [HI Theorem 9.1] we can find v G X® n C°(T) such that 

\\u-v\\ H s {r) <p- 2 ^- s hog^l+p) 

with a = min{A + 1, /x + 1/2} (the integral mean zero condition can be added in a postpro- 
cessing step and does not diminish the convergence rate). This bound holds for s G [0, 1] by 
interpolation. Choosing s = 1/2 and s = l/2 + log _1 (l+p) (where the interpolation norm is uni- 
formly equivalent to the Sobolev-Slobodeckij norm; for details see Section [4] and, in particular, 
Remark 4.1) and combining both bounds proves 

\\u-U hp \\ H}J2{Qh) < (1 + ^/2 log 3/2 (1+p))p -2(a-l/2) log/3(1+p)) 

which is the assertion. □ 
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Remark 3.4. In theory, the additional factor of log 3 / 2 p in the p- estimate by Corollary 3.3 can 



be compensated by choosing v ~ log 3 p. The log 3 ' 2 p perturbation in the error estimate is due to 



the factor (s — 1/2) 3 / 2 in the abstract error estimate of Theorem 3.1 In the case of an open 



smooth s urfa ce (one face), there appears a factor (s — 1/2) l / 2 rather than (s — 1/2) 3//2 ; see 



Theorem ^.7. Therefore, in this case, the p-estimate has only a perturbation of (log p) lf ' 2 which 
can be compensated by choosing v ~ logp instead of v ~ log 3 p as in the case of multiple faces. 
Of course, these are only upper bounds which are not known to be exact and this phenomenon 
is difficult to observe numerically. 

By using /ip-approximation results (for quasi-uniform meshes) from [3] one can easily deduce 
error estimates for the corresponding /ip-version of the discontinuous Galerkin BEM. Neverthe- 
less, in order to keep the presentation simple, we restrict ourselves to the case of arbitrary but 
fixed polynomial degrees. 

Corollary 3.5 (h- version with arbitrary polynomial degree). Let the subspaces X^ p be defined 
by a sequence of meshes with fixed (lowest) polynomial degree p. We assume that the solution u 



of (|2.2|) is of the type (|3.6|) with minimum parameters \i and A and maximum parameters q, s 



and that the regular part u Teg G H 1 ^) H H^ =1 H k (Tj). Then there holds 



\\U - U hp \\ Hl n {Qh) < (1 + V- 112 ) max {^n {fc -l/2 jP+ l/2} ; ^n { A+l/2, M} (J + j logh \f 1+ ^ } 

where 



Pi 

and 



s + 4 + 1/2 if X = pi -1/2, 
s + q otherwise 

\ if p = min{A, 7 - 3}, 
otherwise. 



Proof. By [3, Theorem 7.1] there exists v £ X hp n C°(T) such that 

\ H s (T) < max |/ l mm{fc, P +l}-,^min{A+l, M +l/2}- S ^ + | \ ogh \fl+fo^ . 



U — V 



for < s < min{l, A + 1, fi + 1/2}. Using this result together with the bound from Theorem 3.1 



proves the assertion. □ 



4 Technical results and proof of the main theorem 

Throughout we will make use of the continuity and ellipticity of V proved by Costabel for 
Lipschitz surfaces [6]: 

V : F s_1 (r) -> H S (T), 0<s<l, (4.1) 

(Vv,v) r > \\v\\ 2 H - 1/3(r) Vv G H-^iT). (4.2) 
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We will need several technical results involving fractional order Sobolev norms. So far we have 
defined the Sobolev-Slobodeckij norm. However, standard tool to prove estimates in Sobolev 
spaces is interpolation theory. We therefore have to deal with at least two different definitions 
which, usually, give rise to equivalent norms, cf. Remark 4.1 below. We use the K- method of 
interpolation (see, e.g., [H] for details). For Q C IR n and < s < 1, the following spaces can be 
equivalently defined (i.e., norms are equivalent) 

H s (n) = (l 2 ^),// 1 ^)) , H S {Q) = (L 2 (tt),H^(nf) 
For negative s, we also need spaces defined by interpolation, i.e. 

H%n) = (ii-\n),L\nj) (-i<s<o). 

\ / l+s,2 

Remark 4.1. For fixed domain and fixed order, it is well known that norms in Sobolev spaces 
defined by interpolation and by the double integral (Sobolev-Slobodeckij) are equivalent. In order 
to avoid dependence on a variable domain (elements) we use the equivalence on a fixed domain 
(reference element) by previous transformation. We will need this equivalence for varying order 
s close to 1/2 where it is uniform, cf. Proof of Corollary 4j- 

In order to be transparent, we indicate the type of norm used in the technical lemmas below 
(i.e., specifying "Sobolev-Slobodeckij" or "K-method" or "duality"). 

Lemma 4.2. Let Q be the generic reference element and R a fixed Lipschitz domain (possibly 
Q) with boundary dR. 

(i) (Sobolev-Slobodeckij) For any s G (—1/2, 1/2), and any v G H S (R), there holds 

W v \\h°{r) 5 i/2-\s\ ^ HS{Ry (4 ' 3) 

(ii) (Sobolev-Slobodeckij) For any s G (1/2, 1], and any v G H S {R), there holds 



(4.4) 



IMI^(Q) - h2 ~ 2S W v M Q\\%(Q) ( v G ^(3)' s G [°> X l> K-method), (4.5) 
v\ 2 H s (0 ) ^ h 2 ~ 2s \v o Af Q | 2 - (v €H S (Q), s G [0, 1], Sobolev-Slobodeckij) (4.6) 



\\vh* (dR) < -j===\\v\\ H . w . 

(Hi) For a (bi)linear mapping Mq : Q — > Q onto an element Q G Qh there holds 
\% s{Q) ^h 2 - 2s \\voM Q \\% a ^ 



and 

\\v\\ 2 HS{Q) ~ h 2 - 2s \\vo M Q \\ 2 Hs{0) (v G H s (Q), s G [-1,0], duality). (4.7) 
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Proof. Assertion (i) is shown by |1 1| Lemma 5] and (ii) is Lemma 4.3]. Equivalence (4.6) 
follows from the definition of the semi- norm. For (4.5) and (4.7) in the case of a reference square 
and affine maps see Lemma 2] . The cases of triangular and quadrilateral non-parallelogram 
elements follow similarly. □ 



Lemma 4.3. (K-method) For any s G [— 1, 1] there holds 



QeQ h 



7,H 2 < 

V \\H*(Q) ~ 



v\\ 2 - 



\ V \\H*(T) 

Vu G H S (T 



Vv G H"(T), 



\ r . : /nr,! Vj. 



(4.8) 
(4.9) 



Proof. For complex interpolation these estimates have been analyzed by von Petersdorff |17j ; 
for the K-method see [I], cf. also □ 

Lemma 4.4. (i) 



curiA v 



> 



Q "UH-WiQ) ~ \ V \m/ 2 (Q) 



Vu G # 1/2 (Q) 



(ii) 



curlr : H 1 ' 2 ^) -> (h-V 2 {Y) 
(Hi) (K-method) For s G [0, 1] there holds 



and 



curl rj : H s (Tj 



curl rj : H s (Tj 



,L, 



,L. 



(4.10) 
(4.11) 

(4.12) 
(4.13) 



Proof, (i) is [9l Lemma 4.1] and for (ii) we refer to [3]. In (4.11 ) we simply put a norm for vector 
functions with three components. More precisely, curlr maps onto a vector space of tangential 
fields. For details we refer to [3J. On individual faces, where r, is identified with a subset of 
IR 2 , this tangential space is simply a vector space of functions with two components. In order 
to show (4.12) we recall that Vj is being identified with a subset of IR 2 and, thus, we can extend 
v G L 2 (Tj) by to v° G L 2 (IR 2 ). By Fourier analysis, curl^ : L 2 (IR 2 ) -> (tf-^IR 2 )) 2 , so that 
by the density of smooth functions with compact support in L 2 (Tj) 



curl rj u||#_i (r ) < ||curl H 2 ^H^-in^ < \\v°\\^ m 



Here, the estimate H^H/j--!^-) < II^°IIh- 1 (IR 2 ) ^ or an y ^ e ^o°0-j') with zero extension tp° 
follows from the existence of a bounded extension operator E\ : H 1 (Tj 

sup t: — jj = sup 



MU 2 (r,)- 
o 

H^JR 2 ) 



(4.14) 



H-i(r,-) 



< 



sup 



< 



sup 



H 2 
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Since also curl^ : Hq(Tj) — > (L 2 (Tj)^) , the assertion (4.12) follows by interpolation. The 
boundedness ( 4.13[ ) is proved by interpolation between curlr. : i/ 1 (r j ) — > (L 2 (Tj)) 2 and 
curlrj : L 2 (Tj) — > (-£f -1 (Tj)) 2 . The latter boundedness follows from (4.14) by noting that 

II curl rj v\\ H -i (rj) < || curl rj ulljj-i^). 

□ 

Proposition 4.5 (consistency). Let f G L 2 (F) be given. The DG BEM formulation is consis- 
tent, i.e. the solution u of (2.2) solves 



(«,l)r = 



and 



a h (u,v) = (f,v)i 



Vv G X° hp . 



Proof. By definition u has integral mean zero. It is well known that, for the given geometry and 
/ € L 2 (T), u G ff 1 (r). Therefore u is continuous in the sense of traces and curl/j u = curlr u 
on r, and V curlr u is continuous in the sense of traces, as well. Hence, recalling the signs in 
the definitions of T (3.5), the jump (3.1) and t e ( |3.2[ ), 

a h (u, v) = {V curlr u, curl/j v)r + (Tu, [v]) 7h 

= ^2 ( v curl r u, curl Q v) Q +Y^ ( Tu i M)e 
QeQh ee£ h 

= ^2 {V curl r u,cutI q v)q + ^ (tQ • F curl r u,v) aQ . 



QeQ h 



QeQ h 



Integration by parts on each element proves the assertion, cf. (2.7). 



□ 



Proposition 4.6 (discrete ellipticity) . There exists a constant C > independent of v and h 
such that 



ah(v,v) > (llcurlfculljy-va^) + V^HMHl^)) \M h i/h 



{Qh) VveX° hp . 



Proof. Let v G X? be given. Combining (4.2) and (4.8) we find that there holds 



(V curhi v, curlft v) r > || curl fe 



/2(r) 



> ||curlQu||^_ a 



/2 (Q)' 



(4.15) 
(4.16) 



Transforming Q forth and back onto the reference element Q, denoting the transformed function 
v by v, and employing (4.7), (4.10), (4.6), we obtain 



curlQw|| 2 ¥ _ 1/2(Q) ~ h 3 \\h 1 curlew 



> 



i\*.\2 I I z 
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Together with (4.16) this estimate proves that 



Finally, noting that 



(^curl fc «,curl ft t;) r > ^ \v\ 2 Hl/2{Q) = \v\ 2 Hl/2{Qh y 

QeQ h 

a h (v,v) = {V curl h v,curl h v) r + u([v],[v]) 7h , 



(4.17) 



a combination of (4.15) and (4.17) together with the definition of the norm 
the proof. 

4.1 Proof of the main theorem 



finishes 

□ 



By the discrete ellipticity of the bilinear form a^(-, ■), there exists a unique solution Uh p to (3.3) 
and it remains to prove the error estimate. 

We start with the standard Strang-technique of introducing a discrete function v, in this 
case v E X® p n C°(T), and using the triangle inequality: 



\\u - u hp \\ Hl „ {Qh) < \\u - v\\ Hl/2{Qh) + \\u hp - v\\ Hl/2{QhY 
Since [u] = [v] = on 7^ and {u, l)r = {v, l)r = there holds 

I|U " V hy 2 (Q h ) = \ U ~ V \h^(Qh) ~ ' n " V \h^(F) < \W ~ «llfli/a(r)- 
Here we made use of the bound 

\v\ Hs{ Q h) <C\v\ H s ir) VveH s (T) (a€(0,l]) 



(4.18) 



(4.19) 



(4.20) 



with s = 1/2. For the Sobolev-Slobodeckij semi-norm (employed in this instance to formulate the 
theorem), this follows with C = 1 from its definition. For the norm defined by interpolation one 
uses ( |4.8| ) and an argument from quotient spaces to reproduce the semi-norm on the right-hand 
side. 



To bound the second term on the right-hand side of (4.18) we use the discrete ellipticity of 
the bilinear form a^(-, •), cf. Proposition 4.6 We obtain for any v £ X® p n C°(T) 



n ,| ^ a h (u hp -v,w) 

\\Uhp-v \\ H y2 (0 < sup I, 1, ■ — . 



(4.21) 



any u> G X? 



By Proposition 4.5, the definition of Uh p and the fact that [u] 

^ h P 



on 7/1, there holds for 



ahiuup -v,w) = a h (u - v,w) = (V cur\ h (u - v),curl h w) r + (T(u - v), [w}) lh 

~ ([« - U],rw>7fc + ~ v \> With 

= {VcuY\r{u-v),cur\ h w)r+{T[u-v),[w\) lh . (4.22) 
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The first term on the right-hand side is bounded due to the mapping properties (4.1) of V and 
(Olb of curl r : 



(V curl r (n - v),curl h w) r < ||Fcurl r (u - v) \\ H i/2^\\ curl h w\\ H -i/ 2 ^ 



< 



\u - v 11^1/2^) || curlfe 1011^-1/2^. 



(4.23) 



To bound T{u — v) we employ a trace argument. Let us consider an element Q £ Qh with edge 
e, and a function ip £ H S {Q) with s > 1/2. We map Q onto a fixed reference element Q and 
denote the mapped edge by e. Accordingly we denote the transformed function ip by ip. Then, 
by transformation, using the trace theorem (4.4) and the scaling properties (4.7) with s = and 

h 



(4.6), we find that 



2 

L2(e) 



2 < 
L 2 (e) ~ 



1/2 



/r 



i 2 (Q) 



+ ^ s "Vl 2 



H°(Q) 



Applying this bound to the vector case, and recalling the definition (3.5) of T, yields (with 
Sobolev-Slobodeckij semi-norm) 

1 

I hi II 1 / mini I in « i \ 1 1 

(Q) 



< 



li 2 (e) ~ s _ 1/2 

and, summing over e € £h, 



h- 1 1| 7 curlr (u - v)f mQ) + ^ 2s ~ 1 |F curl r ( 



u — V) 



\T(u — v) 



< 



1 



1/2 



h- 1 1| y curlr {u - u)|||a (r) + /i 2s_1 |y curl r 



■u — f 



By the mapping properties (4.1) of V, the norm estimate (4.9) and the boundedness (4.12) of 
curlr- (j = 1, • ■ ■ j L) we have that 



\V curl T (u - v)\\ 2 L 2 {r) < ||curl r (u-f)||^_ 1(r) < || curl rj (u - v) 

L 



Ifl-l(Ii) 



< 



u — V 



\u — v\ 



L2(r)- 



(4.25) 



Similarly, first employing (4.20), then the mapping properties (4.1) of V, the decomposition of 
norms (4.9), the norm equivalence (4.3), the boundedness of the surface curl (4.13), and (4.8), 
we obtain 



\V curl r (n - v)\ 2 Hs{Qh) < |Fcurl r (u - v)\ 2 H3{r) < || curl r (u - u)||^._i 



(T) 



< II curl rj ( 



u — v) 



< 



1 



1 y 

(s - 1/2)2 



u ~~ u lllr»(r 3 -) ~ 



[s - 1/2)2 Z_. 

1 



Y, II curl rj ( 



(5-1/2)2 



u — V 



(4.26) 
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Here, we used implicitly that the Sobolev-Slobodeckij norm (started with) is uniformly equivalent 
to the interpolation norm for s close to 1/2 so that the continuity (4.13) of curlr with respect 
to the interpolation norm is applicable. In the final term we can switch back to the Sobolev- 
Slobodeckij norm. 



Combining (4.24), (4.25) and (4.26), we obtain 



1 



\\T(u-v)\\ L 2 (7h) 



< 



1 



" 1/2) 



= (V 1/2 I 

1/2) V 



u-v\\ L 2 {T) + 



s - 1/2 



(4.27) 



Combination of pT^, (|422|), (|423|) and f07|) proves 
11?%, 



< 



|u - u|| H i/2 (r ) + 



1 



vM* - V2) 



[h- l/2 \\u - v\\ L 2 {T) + -— — ||u - v\\ H s {v) 



Hv (Qh) 

for any u G n C°(r). By ( |4.18[ ) and ( |4.19[ ) this concludes the proof. 
4.2 The problem on an open surface 

In the case of an open polyhedral surface the previous analysis carries over with few changes. 
The energy space (where one looks for the solution u of (2.1)) then is H l l 2 (T). Moreover, the 
hypersingular and weakly singular operators W and V are elliptic on H l l 2 (T) and i/ -1 / 2 (r), 
respectively, and the continuity (4.1) becomes 



V : H a - L (T) -> H S (T), < s < 1. 



(4.28) 



Since in this case the hypersingular operator is invertible one does not ask for u to have integral 
mean value zero. Instead, its trace on the boundary dT of T vanishes. Accordingly, functions 
of the discrete space Xh p must comply with this homogeneous boundary condition. In the 
corresponding discrete bilinear form a,h(-,-), jumps on edges lying on dT become traces. 

To be precise, let us recall the notation Eh for the set of edges of Qh- We distinguish between 
edges E\ which are subsets of the boundary dT, and the rest E h , 



£h = U £ % h . 



For e G Si with e = Q 1 n Q 2 , QuQz G Q h , 



as before, and for e G 



[t>J| e := vie- 
Accordingly we define unit tangential vectors on edges, 



C Qile 



ifeG^, 
ifeG^ 



(4.29) 
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with t being the unit tangential vector along dT. The discontinuous Galerkin boundary element 
method then reads: Find u^p G X^ p such that 

ah(u hp , v) = (f, v) r W G X hp , 



and, instead of Theorem 3.1, there holds the error estimate 
\\u-U hp \\ H y 2{Qh) <C 

il — v\\h1/2(t) + 



inf f I 
rim L 



1 



vex hp nH%(F) 



/l/2( 



1/2) 3 / 2 



(h- 1/2 \ 



u — V 



Il2 ( d + 



h s - l ' 2 \ 



u — V 



H s (r) 



(4.30) 



with s G (l/2,min{r, 1 — 5}] for r < 1 (since u .H^r) in this case) and norm 
Let us note the few changes which are necessary in the proofs. 



1. Proposition 4.5 In the case of an open surface we have u G H S (T) for any s < 1 so that 
the jumps of u on 7^ including its trace on dT vanish, and V curlr u is also continuous. 
The rest of proving consistency is identical. 



2. Proposition 4.6 In (4.15) one replaces the H- l / 2 (T) with the H' 1 / 2 ^) -norm and notes 



that the injection of H 1 I 2 (T) in H 1//2 (r) is continuous. 



3. Proof of (4.30). The steps from Section 4.1 carry over by noting that [u] = u = on dT 
and [v] = v = on dT for any v G Xh p , and by replacing some "non-tilde" norms with 
their "tilde" counterparts. 

Theorem 4.7. In the case of a single face the error estimate can be improved to 



\u-U hp \\ H yi {Qh) 



< C 



inf ( 

veX hp nHi(T) L 



u — v\ 



J H" 1 /2(r) 



+ -T^rm 1/2||n " vh2{r) + ^ 1/2||n " vhs ^y 



\M« - 1/2) 

Proof. This estimate follows as discussed above and by replacing (4.26) with the estimate 

\h'(t) ~ II curl r(^ - v)\\%s-i (r) Z " ■ ••" 2 



\V curl r (u-v)\ 2 HS(Qh) < |^curl r ( 



u — V ) 



u — v\ 



H s (r) 



where we used the continuity (4.28) of V and the continuity (4.12) of curlr on the single face 
r. Then, instead of (4.27), we obtain 

M*- V2) 



irn»-'-)|i^,-, ; ,,< 1/2) (ft~ 1/a i 



u — V 



\l\v) + 



it — u 



H-(r);- 



and the result follows. 



□ 
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Remark 4.8. The improved estimate from Theorem \4- 7| leads to the same convergence order 



of the h-version as stated by Corollary 3.2 (though r < 1 by the reduced regularity). For the 
p-version, however, it improves the estimate by a factor of log _1 (l + p), cf. Corollary 



3.3 and 



Remark 3.4 



5 Numerical results 



We consider the model problem (2.1) with / = 1 on V = (0, 1) x (0, 1) and use uniform meshes 
Qh on T consisting of squares. 

Since the exact solution u of ( |2.1[ ) is unknown, the error 

W ~ u hp \ 2 Hl/2(Qh) + u\\[u hp \\\ 2 L 2 ilh) 

cannot be computed directly, except for || [uh p ] Hl 2 ^) which is straightforward to implement. As 
in previous papers on non-conforming approximations of hypersingular operators (see, e.g., [5]) 
we approximate an upper bound to the semi- norm \u — u hp\H 1 / 2 (Q h )- First, note that there holds 

\u - u hp \ 2 Hl/2 , Q , < (V curbi(u - u hp ), cur\ h (u - u hp ))r, 

cf. ( |4l7| . We find that 

\u - u hp \ 2 H i/2 ( Qfe) < {V curl r u, curl r u) r + (V curU u hp , curl^ u hp ) r - 2{V curl r u, curl^ u hp ) v . 

Let us analyze the three terms on the right-hand side, 
(i) Since u G H l l 2 (Y), there holds 

(V curlp u, curlr u)r = (Wu,u)r- 

The energy norm (squared) (Wu,u)y of u can be approximated through extrapolation, in the 
following denoted by ||u||e X > see |8j. 



(ii) By construction of Uh p (3.3) 



(V curlft u hp , curl^ u hp ) T = (f, u hp ) T - v\\ [u hp ] \\ 2 L2 



(iii) By consistency (see Proposition 4.5), 



(V curlr u, cur l h u hp ) T = (/, u hp ) r - {Tu, [u hp ])~ jh . 
For our example, u 6 ^ 1_e (r) for any e > 0, so that Tu G H 1 ' 2 ~ t ( , yt l ) for any e > 0. By duality, 

(Tu, [u hp ]) lh < \\Tu\\ m/ 2-, {lh) \\[u hp \\\ H -y2+, {lh) . 
Combining (i)-(iii) we conclude that for any e > there holds (with a constant depending on e) 
\u-u hp \ 2 Hl/2{Qh) < \ {Wu,u) T - (f,u hp ) r \ +v\\ [uh P ]\\ 2 L 2 {jh) + \\[uh P }\\ H -i/2+e {yh) . 
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Since u 6 H 1 / 2 € {^h) (and considering the singularities of u) we expect that 



MlH-V 2+ e( 7 ,) < ^- 2 V 2(1 - 2£) 



and 



are of the same order (note the different exponents of the norms), so that for our numerical 
experiments we take 



MIL - (f> u h P )r\ 1/2 + \\[uhp]\\L2( Jh ))/\\u\ 



as a computable and reasonable measure for an upper bound of the error \\u 
normalized by ||«|| jji/2fln ~ IMUx- 

Below we present numerical results for the two contributions 



u hp\\H l / 2 (Q h ) 



U 



(f, u hp )r\ 1/2 /\\u\ 



(referred to as "i^ 1 / 2 " error in the figures) and 



[Uhpl 



1 1/2 

Il2 (7) 



/[Ml 



(referred to as "L error) which measures the jumps, a measure for non-conformity. 



We first study the influence of v. In Figure 5.1 four different approximations are presented, 
all on the mesh of 25 elements and polynomial degree p = 3. The first three pictures (a)- 
(c) show Uh p for increasing v (y = 0.1, 1, 10) and the last picture (d) shows the conforming 
approximation. One sees that with increasing v discontinuities quickly disappear, and Uh p with 
v = 10 is visually not distinguishable from the conforming counterpart. Note, in particular, 
how the trace of Uh p on the boundary of T approaches when v increases. Nevertheless, even 
for v = 100, discontinuities are still present and converge at the predicted rate, see Figure 5.7 
below. 

Now let us study convergence rates. In our case of an open surface, the strongest singularities 
in ( 3.6 ) are edge singularities with fi = 1/2 so that the h- version for any polynomial degree should 

see [21 [3], and HU for numerical results. 



converge like h 1 ' 2 and the jo-v ersion like p 



In Figures 5.2, 5.4 and 5.6 we present on double-logarithmic scales the relative "iJ 1 / 2 " errors 
with v = 1,20,100, respectively, for the polynomial degrees p = 1,3,5 in all cases. We also 
show the errors for the conforming method (p = 1) and, for comparison, the curve h 1 / 2 . The 



latter curve indicates the expected rate of convergence (cf. Corollaries 3.2 and 3.5) and, indeed, 



the numerical results show at least this convergence rate. Note, however, that in Figure 5.2 
the curve for p = 1 is below the ones for p = 3,5, whereas for a conforming method increasing 
p reduces the error. In this case v = 1 is small and we explain this phenomenon by the fact 
that we are dealing with a saddle point formulation where on cannot expect convergence like 
that of a projection method (as it is the case with conforming approximations of this operator). 
Indeed, when v increases, the expected order of the curves is re-established. First, for v = 20 
in Figure 5.4 we have a pre-asymptotic super-convergence and then, for v = 100 in Figure [5T6[ 
the curves for the different degrees have the right order (highest degree gives smallest error). 
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In Figures 5J3 5J5 
of the right 
In Figures 



and 



5.7 the corresponding "L errors of the jumps are given. They are 



all of the right con verg ence order h 1 ' 2 and in the right order. 

and 5.9 we plot the errors ("ff 1 / 2 " resp. "L 2 ") obtained by the p- version 
and for different values of v [y = 1, 10, 50, 100) along with the curve p~ l which represents the 



expected convergence rate, cf. Corollary |3.3| For the given range of degrees of freedom, smaller 
v does not result in the expected convergence order of the V 2 " errors (at least not of our 
upper bound) and we conclude that we are in a pre-asymptotic range. For higher u [y = 50, 100) 
the optimal order is observed. In contrast, the "L 2 " convergence of the jumps is the optimal 
one for all v and, not surprisingly, larger v leads to smaller errors. 
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(c) v = 10 (d) conforming 



Figure 5.1: Approximation with 5x5 elements and p = 3. 
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1 




DG h-BEM (nu=1, p=1), H A {1/2}-error — +— 

DG h-BEM (nu=1, p=3), IT{1/2}-error + 

DG h-BEM (nu=1, p=5), H A {1/2}-error — +- 



100 1000 10000 

dof 

Figure 5.2: /i-version with degrees p = 1,3,5 and v = 1, relative "^V 2 " errors. Comparison 
with conforming BEM. 



1 



<D 0.1 



0.01 




h A {1/2} 

DG h-BEM (nu=1, p=1), L2-error(jump) 
DG h-BEM (nu=1, p=3), L2-error(jump) 
DG h-BEM (nu=1, p=5), L2-error(jump) 



100 



1000 



dof 



10000 



Figure 5.3: /i-version with degrees p = 1,3,5 and v = 1, relative "L 2 " errors 



20 



h A {1/2} 

conforming h-BEM (p=1), error in energy norm — i — 
DG h-BEM (nu=20, p=1), H A {1/2}-error --+- - 

DG h-BEM (nu=20, p=3), H A {1/2}-error + 

DG h-BEM (nu=20, p=5), H A {1/2}-error — +- 




1 



0.1 



0.01 



0.001 



h A {1/2} 

DG h-BEM (nu=20, p=1), L2-error(jump) — +- 

DG h-BEM (nu=20, p=3), L2-error(jump) +■■■ 

DG h-BEM (nu=20, p=5), L2-error(jump) ■■■+■• 



100 



1000 



dof 



Figure 5.5: /i-version with degrees p = 1,3,5 and v = 20, relative "L 2 " errors 



21 



h A {1/2} 

conforming h-BEM (p=1), error in energy norm — i — 

DG h-BEM (nu=100, p=1), H A {1/2}-error --+- - 

DG h-BEM (nu=100, p=3), H A {1/2}-error + 

DG h-BEM (nu=100, p=5), IT{1/2}-error — +- 




100 



1000 



10000 



dof 



Figure 5.6: /i-version with degrees p = 1,3,5 and v = 100, relative "iJ 1 / 2 " errors. Comparison 
with conforming BEM. 



<D 0.01 



0.001 



0.0001 



1^ 

+ , ^—-H 

+.... +■ + 



h A {1/2} 

DG h-BEM (nu=100, p=1), L2-error(jump) — +— 

DG h-BEM (nu=100, p=3), L2-error(jump) + 

DG h-BEM (nu=100, p=5), L2-error(jump) ■■-+-■ 



100 



1000 



10000 



dof 



Figure 5.7: /i-version with degrees p = 1, 3, 5 and v = 100, relative "L 2 " errors. 
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10 100 

dof 



Figure 5.8: p-version with 4 elements and v = 1, 10,50, 100, relative "if 1 / 2 " errors. 
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P A {-1} 






DG p-BEM (nu=1), L2-error(jump) — +— 






DG p-BEM (nu=10), L2-error(jump) + 






DG p-BEM (nu=50), L2-error(jump) -■-+■-■ 






DG p-BEM (nu=100), L2-error(jump) — ■+■— 





10 100 

dof 

Figure 5.9: p-version with 4 elements and v = 1, 10, 50, 100, relative "L 2 " errors. 
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